Photonic crystal ribbon-beam traveling wave amplifier

ABSTRACT

A RF amplifier includes a RF input section for receiving a RF input signal. At least one single-sided slow-wave structure is associated with the RF interaction section. An electron ribbon beam that interacts with the RF input supported by the at least one single-sided slow-wave structure so that the kinetic energy of the electron beam is transferred to the RF fields of the RF input signal, thus amplifying the RF input signal. A RF output section outputs the amplified RF input signal.

PRIORITY INFORMATION

This application claims priority from provisional application Ser. No. 60/483,852 filed Jun. 30, 2003, which is incorporated herein by reference in its entirety.

This invention was made with government support under Grant No. F49620-03-1-0230 awarded by the Air Force. The government has certain rights in the invention.

BACKGROUND OF THE INVENTION

The invention relates to the field of optical communication, and in particular to a photonic crystal ribbon-beam traveling wave amplifier.

The third-generation (3G) wireless communication standards call for hardware-based upgrade to the second-generation (2G) Global Systems for Mobile Communications using Wideband Code Division Multiple Access (W-CDMA) and Universal Mobile Telephone System (UMTS) as well as software-based upgrade to 2G Code Division Multiple Access (CDMA). The 3G wireless communications require amplifiers operating frequencies that are 1.12 to 3 times that of present frequencies, which are in the range from 900 MHz to 1700 MHz.

In general, the bandwidth of a transmitter, which is the most important figure of merit, increases with the central frequency of the amplifier. However, the number of transmitting towers must increase as the square of the central frequency, while keeping the power of the transmitting tower at a constant. This is because the distance between two adjacent transmitting towers is inversely proportional to the frequency. For example, if 1-GHz transmitting towers have a spacing of 10 miles, then 2-GHz transmitting towers must have a spacing of 5 miles. In other words, four 2-GHz transmitting towers are required to cover 100 square miles, whereas only one 1-GHz transmitting tower is needed for the same area.

Moreover, the total RF power per unit area increases with increasing data rate. For example, 3G wireless networks are expected to have considerably higher data rate than 2G wireless networks. As a result, the number of power amplifiers increases more dramatically than the square of the carrier frequency.

In order for the third-generation (3G) and future wireless communications to be a viable business, it is essential for the telecommunication equipment industry to provide ultra-low-cost amplifiers.

At present, most wireless base stations are powered by solid-state power amplifiers, which operate with efficiencies in the 8-12% range. The cost of solid-state amplifier is about $100/Watt. The cost of power amplifiers per base station at 1.5 kW is $150,000. For the RF power part of a wireless base station, the operating cost is comparable to the capital cost, because of the low operating efficiencies and heat removal.

Conventional helix traveling wave tubes (TWTs), which are not employed in any existing wireless base stations, cannot meet the ultra-low-cost requirement set by any potential third-generation wireless infrastructure provider. For example, a 100 W, 2 GHz conventional helix TWT costs $20K a piece or more.

There is a need to develop high efficiency, low-cost microwave amplifiers for 3G and future wireless base stations.

SUMMARY OF THE INVENTION

According to one aspect of the invention, there is provided a RF amplifier. The RF amplifier includes a RF input section for receiving a RF input signal. At least one slow-wave structure associated with the RF interaction section. An electron ribbon beam interacts with the RF input supported by the at least one slow-wave structure so that the kinetic energy of the electron beam is transferred to the RF fields of the RF input signal, thus amplifying the RF input signal. A RF output section outputs the amplified RF input signal.

According to another aspect of the invention, there is provided a method of forming a RF amplifier. The method includes forming RF input section for receiving a RF input signal. Also, the method includes forming al least one photonic crystal for operational control if necessary. An electron ribbon beam is formed that interacts with the RF input supported by the at least one slow-wave structure so that the kinetic energy of the electron beam is transferred to the RF fields of the RF input signal, thus amplifying the RF input signal. Furthermore, the method includes forming a RF output section that outputs the amplified RF input signal.

BRIEF DESCRIPTION OF THE DRAWINGS

FIGS. 1A-1B are schematic diagrams of a photonic crystal ribbon-beam traveling wave amplifier (PCRB TWA) in accordance with the invention;

FIG. 2 is a schematic diagram illustrating a detailed depiction of a single-sided slow-wave structure;

FIG. 3 is a schematic diagram illustrating a detailed depiction of a double-sided slow-wave structure;

FIGS. 4A-4B show plots of the dispersion functions D_(a)(ω,k) for the anti-symmetric modes as they vary with frequency at several values of phase shifts;

FIGS. 5A-5B show plots of the dispersion functions D_(S)(ω,k) for the symmetric modes as they vary with frequency at several values of phase shifts;

FIG. 6 is a graph illustrating eigenfrequencies for the two lowest bands of anti-symmetric modes and the two lowest bands of symmetric modes;

FIG. 7 is a graph illustrating the comparison between the eigenfrequencies for the two lowest anti-symmetric modes with the same system parameters as in FIG. 6;

FIG. 8 is a schematic diagram illustrating a ribbon electron beam interacting with a single-sided PC slow-wave structure;

FIGS. 9A-9B are graphs illustrating the temporal linear growth rate of the lowest anti-symmetric mode as a function of the wavenumber;

FIGS. 10A-10B are graphs illustrating the spatial linear growth rate of the lowest anti-symmetric mode as a function of angular frequency;

FIG. 11 is a graph demonstrating the maximum temporal growth of the lowest anti-symmetric mode as a function of normalized beam current;

FIG. 12 is a graph demonstrating the maximum spatial growth of the lowest anti-symmetric mode as a function of normalized beam current;

FIGS. 13A-13B is a graph demonstrating the dispersion characteristics in a two-dimensional dielectric photonic crystal;

FIG. 14 is a schematic diagram of the design of the PCRB TWA for 3G wireless base stations;

FIG. 15 is a graph of the eigenfrequency versus the phase shift for the two lowest bands in the 200 W, 1950 MHz, 3% bandwidth PCRB TWA structure;

FIG. 16 is a graph of the normalized real wavenumbers Rek_(z)L versus the normalized frequency for the operating mode in the 200 W, 1950 MHz, 3% bandwidth PCRB TWA structure;

FIG. 17 is a graph of the normalized imaginary wavenumbers |Imk_(z)L| versus the normalized frequency for the operating mode in the 200 W, 1950 MHz, 3% bandwidth PCRB TWA structure;

FIG. 18 is a graph of the real phase shift Rek_(z)L versus the frequency f for the operating mode;

FIG. 19 is a graph of the gain G versus the frequency f for the operating mode; and

FIG. 20 is a graph demonstration the dispersion characteristics for 4 lowest-order TM modes in a photonic crystal.

DETAILED DESCRIPTION OF THE INVENTION

The invention is a novel amplifier that employs two emerging technologies, namely, photonic crystals and low-density ribbon electron beams, in otherwise a conventional vacuum tube millimeter wave amplifier.

FIGS. 1A-1B shows a schematic diagram of a photonic crystal ribbon-beam traveling wave amplifier (PCRB TWA) 2 in accordance with the invention. FIG. 1A shows a double sided PCRB TWA 2 that includes a ribbon electron beam 10 propagating in the z-direction from the emitter 14 and extending out, wiggler magnets 8 for beam focusing, a photonic crystal (PC) slow-wave structure 12 with metallic or dielectric rods and plates, and RF input 4 and output 6 sections. As the electron beam 10 interacts with the RF input 4 supported by slow-wave structure 12, the kinetic energy of the electron beam is transferred to the RF fields, amplifying the RF signal 16. The amplified RF signal 18 exits the amplifier at the RF output 6, and the spent electron beam is collected down stream.

FIG. 1B shows a single sided PCRB TWA 20 that includes a ribbon electron beam 28 propagating in the z-direction from the emitter 32 and extending out, wiggler magnets 26 for beam focusing, a photonic crystal (PC) slow-wave structure 30 with metallic or dielectric rods and plates, and RF input 22 and output 24 sections. As the electron beam 28 interacts with the RF input signal 34 supported by slow-wave structure 30, the kinetic energy of the electron beam 28 is transferred to the RF fields, amplifying the RF signal 34 to produce the amplified signal 36.

Unlike the round beams in conventional helix or coupled-cavity TWTs, the ribbon electron beam in the PCRB TWA will reduce the magnetic field required for beam focusing, reduce the loading in the amplifier, increase the amplifier efficiency, and improve the amplifier linearity and bandwidth.

The ribbon beam will have an aspect ratio of 1 to 10, which effectively lowers the beam perveance (or space-charge) by a factor of 10 in comparison with a round beam. Based on the well-known empirical scaling law, the efficiency of the PCRB TWA is expected to be as high as 80%. Furthermore, because the effective beam perveance is small, the interaction between the electrons and slow-wave structure is expected to be high. Consequently, a high degree of linearity of the amplifier is expected in the high-efficiency operation.

Consider a monochromatic electromagnetic wave propagating in the two-dimensional single-sided slow-wave structure 38 as shown in FIG. 2. The slow-wave structure 38 consists of a metal plate at x=0 and corrugated vanes 40 located between x=b and d. The period of corrugation is L, and the width of each vane is a. For simplicity, either the periodic metal or dielectric rods that are placed at x≧d in the PCRB TWA or the variations in the y-direction is considered. Moreover, the transverse magnetic (TM) modes are of importance with the field distributions E(x,z,t)=e ^(−iαx) [E _(x)(x,z)ê _(x) +E _(z)(x,z)ê _(z)],  Eq. 1 B(x,z,t)=e ^(−iαx) B _(y)(x,z)ê _(y).  Eq. 2 where E(x,z,t) is the electric field, E_(x)(x,z) is the electric in the x-direction, E_(z)(x,z) is the electric field in the y-direction, B(x,z,t) is the magnetic field, and treat B_(y)(x,z) as the generating function. The wave equation can be expressed as

$\begin{matrix} {{{\left( {\frac{\partial^{2}}{\partial x^{2}} + \frac{\partial^{2}}{\partial z^{2}}} \right){B_{y}\left( {x,z} \right)}} + {\frac{\omega^{2}}{c^{2}}{B_{y}\left( {x,z} \right)}}} = 0.} & {{Eq}.\mspace{14mu} 3} \end{matrix}$ where ω is the frequency of the wave and c is the speed of the light. In cgs units, the electric field can be expressed as

$\begin{matrix} {{{E_{x}\left( {x,z} \right)} = {\frac{c}{i\;\omega}\frac{\partial}{\partial z}{B_{y}\left( {x,z} \right)}}},} & {{Eq}.\mspace{14mu} 4} \\ {{E_{z}\left( {x,z} \right)} = {{- \frac{c}{i\;\omega}}\frac{\partial}{\partial x}{{B_{y}\left( {x,z} \right)}.}}} & {{Eq}.\mspace{14mu} 5} \end{matrix}$

Expressing the generating function as a Bloch-wave function of the form

$\begin{matrix} {{{B_{y}\left( {x,z,t} \right)} = {{\mathbb{e}}^{i{({{k_{z}z} - {\omega\; t}})}}{\sum\limits_{n = {- \infty}}^{\infty}{{u_{n}(x)}e^{{i2}\;\pi\;{{nz}/L}}}}}},} & {{Eq}.\mspace{14mu} 6} \end{matrix}$ one can rewrite the wave equation as

$\begin{matrix} {{{\sum\limits_{n = {- \infty}}^{\infty}{\left\{ {\frac{\partial^{2}u_{n}}{\partial x^{2}} + {\left\lbrack {\frac{\omega^{2}}{c^{2}} - \left( {k_{z} + \frac{2\;\pi\; n}{L}} \right)^{2}} \right\rbrack u_{n}}} \right\}{\mathbb{e}}^{{i2}\;\pi\;{{nz}/L}}}} = 0},} & {{Eq}.\mspace{14mu} 7} \\ {or} & \; \\ {{{\sum\limits_{n = {- \infty}}^{\infty}{\left\{ {\frac{\partial^{2}u_{n}}{\partial x^{2}} + {P_{n}^{2}u_{n}}} \right){\mathbb{e}}^{{i2}\;\pi\;{{nz}/L}}}} = 0},} & {{Eq}.\mspace{14mu} 8} \end{matrix}$ where u_(n) is the Floquet amplitude, and p_(n) is an effective wavenumber.

$\begin{matrix} {p_{n}^{2} = {{p_{n}^{2}\left( {\omega,k_{z}} \right)} = {\frac{\omega^{2}}{c^{2}} - {\left( {k_{z} + \frac{2\;\pi\; n}{L}} \right)^{2}.}}}} & {{Eq}.\mspace{14mu} 9} \end{matrix}$ Note that for

$\begin{matrix} {{{\frac{\omega^{2}}{c^{2}} - \left( {k_{z} + \frac{2\;\pi\; n}{L}} \right)^{2}} < 0},} & {{Eq}.\mspace{14mu} 10} \end{matrix}$ p_(n) is imaginary.

In the corrugated-vane region, the usual approximation is adopted

$\begin{matrix} {{{B_{y}\left( {x,z,t} \right)} \cong {{iE}_{0}{\mathbb{e}}^{i{({{k_{z}{sL}} - {\omega\; t}})}}{\cos\left\lbrack {\frac{\omega}{c}\left( {x - d} \right)} \right\rbrack}}},} & {{Eq}.\mspace{14mu} 11} \end{matrix}$ E _(x)(x,z,t)≅0,  Eq. 12

$\begin{matrix} {{E_{z}\left( {x,z,t} \right)} \cong {E_{0}{\mathbb{e}}^{i{({{k_{z}{sL}} - {\omega\; t}})}}{\sin\left\lbrack {\frac{\omega}{c}\left( {x - d} \right)} \right\rbrack}}} & {{Eq}.\mspace{14mu} 13} \end{matrix}$ for |z−sL|<a/2 and b<x<d. Here, s=0, ±1, ±2, . . . Note that Eq. 13 assures E_(z)|_(x=d)=0.

In the drift region with 0<x<b, Eq. 8 implies

$\begin{matrix} {{{\frac{\mathbb{d}^{2}u_{n}}{\mathbb{d}x^{2}} + {p_{n}^{2}u_{n}}} = 0},} & {{Eq}.\mspace{14mu} 14} \end{matrix}$ which has a general solution of the form u _(n)(x)=A _(n) sin(p _(n) x)+B _(n) cos(p _(n) x),   Eq. 15 where A_(n) and B_(n) are constants. Because E_(z)|_(x=0)=0, one must have A_(n)=0,  Eq. 16

or u _(n)(x)=B _(n) cos(p _(n) x).  Eq. 17

Therefore, the RF fields in the drift region can be expressed as

$\begin{matrix} {{{B_{y}\left( {x,z,t} \right)} = {{\mathbb{e}}^{i{({{k_{z}z} - {\omega\; t}})}}{\sum\limits_{n = {- \infty}}^{\infty}{B_{n}{\cos\left( {p_{n}x} \right)}e^{{i2}\;\pi\;{{nz}/L}}}}}},} & {{Eq}.\mspace{14mu} 18} \\ {{{E_{x}\left( {x,z,t} \right)} = {{\mathbb{e}}^{i{({{k_{z}z} - {\omega\; t}})}}{\sum\limits_{n = {- \infty}}^{\infty}{\frac{c}{\omega}\left( {k_{z} + \frac{2\;\pi\; n}{L}} \right)B_{n}\mspace{11mu}{\cos\left( {p_{n}x} \right)}e^{{i2}\;\pi\;{{nz}/L}}}}}},} & {{Eq}.\mspace{14mu} 19} \\ {{E_{z}\left( {x,z,t} \right)} = {{\mathbb{e}}^{i{({{k_{z}z} - {\omega\; t}})}}{\sum\limits_{n = {- \infty}}^{\infty}{\frac{c}{i\;\omega}p_{n}B_{n}\mspace{11mu}{\sin\left( {p_{n}x} \right)}e^{{i2}\;\pi\;{{nz}/L}}}}}} & {{Eq}.\mspace{14mu} 20} \end{matrix}$ for 0<x<b. The amplitudes B_(n) in the drift region are related to the amplitude E₀ in the corrugated-vane region by the continuity of the electric field at x=b, i.e.,

$\begin{matrix} {{{\mathbb{e}}^{{{\mathbb{i}}\; k},z}{\sum\limits_{n = {- \infty}}^{\infty}\;{\frac{c}{i\;\omega}p_{n}B_{n}{\sin\left( {p_{n}b} \right)}{\mathbb{e}}^{{\mathbb{i}2\pi}\;{{nz}/L}}}}} = \left\{ \begin{matrix} {{E_{0}{\mathbb{e}}^{{\mathbb{i}}{({ksL})}}{\sin\left\lbrack {\frac{\omega}{c}\left( {b - d} \right)} \right\rbrack}},} & {{{{if}\mspace{14mu}{{z - {sL}}}} < {a/2}},} \\ {0,} & {{otherwise},} \end{matrix} \right.} & {{Eq}.\mspace{14mu} 21} \end{matrix}$ where s=0, ±1, ±2, . . . Solving Eq. (3.1.21) for B_(n) yields

$\begin{matrix} {B_{n} = {\frac{\omega}{ic}\frac{{abE}_{0}}{L}{\frac{{\sin\left\lbrack {\frac{\omega}{c}\left( {d - b} \right)} \right\rbrack}{\sin\left\lbrack {\left( {k_{z} + \frac{2\pi\; n}{L}} \right)\frac{a}{2}} \right\rbrack}}{\left( {k_{z} + \frac{2\pi\; n}{L}} \right)\frac{a}{2}}.}}} & {{Eq}.\mspace{14mu} 22} \end{matrix}$

The vacuum dispersion relation can be derived with matching the Poynting flux at x=b. For present purposes, it is useful to introduce the admittance defined by

$\begin{matrix} {{A^{\pm}\left( {\omega,k_{z}} \right)} \equiv {\frac{{{\int_{{- L}/2}^{L/2}{E_{z}^{*}B_{y}\ {\mathbb{d}z}}}}_{x = {b \pm 0}}^{\;}}{{{\int_{{- L}/2}^{L/2}{E_{z}{\mathbb{e}}^{{- {\mathbb{i}}}\; k_{z}z}\ {\mathbb{d}z}}}}_{x = {b \pm 0}}^{2}}.}} & {{Eq}.\mspace{14mu} 23} \end{matrix}$ Substituting Eqs. 11 and 13 to Eq. 22, one obtains

$\begin{matrix} {{A^{+}\left( {\omega,k_{z}} \right)} = {\frac{1}{ia}{\cot\left\lbrack {\frac{\omega}{c}\left( {d - b} \right)} \right\rbrack}{\frac{\left( \frac{k_{z}a}{2} \right)^{2}}{\sin^{2}\left( \frac{k_{z}a}{2} \right)}.}}} & {{Eq}.\mspace{14mu} 24} \end{matrix}$ Similarly, substituting Eqs. 18, 20 and 22 into Eq. 23, one obtains

$\begin{matrix} {{A^{-}\left( {\omega,k_{z}} \right)} = {\sum\limits_{n = {- \infty}}^{\infty}\;{\frac{i\;\omega}{{cLp}_{n}}{\cot\left( {p_{n}b} \right)}\frac{\left( \frac{k_{z}a}{2} \right)^{2}}{\sin^{2}\left( \frac{k_{z}a}{2} \right)}{\frac{\sin^{2}\left\lbrack {\left( {k_{z} + \frac{2\pi\; n}{L}} \right)\frac{a}{2}} \right.}{\left\lbrack {\left( {k_{z} + \frac{2\pi\; n}{L}} \right)\frac{a}{2}} \right\rbrack^{2}}.}}}} & {{{Eq}.\mspace{14mu} 25}\;} \end{matrix}$ By setting A ⁺(ω,k _(z))=A ⁻(ω,k _(z)),  Eq. 26 one arrives at the vacuum dispersion relation

$\begin{matrix} \begin{matrix} {{D_{a}\left( {\omega,k_{z}} \right)} \equiv {{\frac{\omega\;{ab}}{cL}{\sum\limits_{n = {- \infty}}^{\infty}\;{\frac{\cot\left( {p_{n}b} \right)}{p_{n}b}\frac{\sin^{2}\left\lbrack {\left( {k_{z} + \frac{2\pi\; n}{L}} \right)\frac{a}{2}} \right\rbrack}{\left\lbrack {\left( {k_{z} + \frac{2\pi\; n}{L}} \right)\frac{a}{2}} \right\rbrack^{2}}}}} +}} \\ {{\cot\left\lbrack {\frac{\omega}{c}\left( {d - b} \right)} \right\rbrack} = 0} \end{matrix} & {{Eq}.\mspace{14mu} 27} \end{matrix}$ for the electromagnetic wave in the single-side slow-wave structure. As pointed out earlier, when inequality of Eq. 10 holds, p_(n) is imaginary and

$\begin{matrix} {\frac{\cot\left( {p_{n}b} \right)}{p_{n}b} = {- {\frac{\coth\left( {{p_{n}}b} \right)}{{p_{n}}b}.}}} & {{Eq}.\mspace{14mu} 28} \end{matrix}$

In general, Eq. 27 must be solved numerically.

As an alternative to the single-side slow-wave structure, a monochromatic wave propagation in a double-side slow-wave structure 42 is considered, as shown in FIG. 3. The bottom of the double-side symmetric slow-wave structure is a mirror image of the top of the single-side slow-wave structure shown in FIG. 2. In the double-side slow-wave structure, however, there is no metal at x=0.

In the double-side slow-wave structure, the vacuum dispersion relation for the anti-symmetric wave propagation with E(−x,z,t)=−E(x,z,t) is the same as the one given in Eq. 27.

However, the vacuum dispersion relation for the symmetric wave propagation with E(−x,z,t)=E(x,z,t) has a difference expression. Paralleling the analysis in Section 3.1, and taking A_(n)≠0 and B_(n)=0 in Eq. 15, one can show that it is given by

$\begin{matrix} \begin{matrix} {{D_{s}\left( {\omega,k_{z}} \right)} = {{{- \frac{\omega\;{ab}}{cL}}{\sum\limits_{n = {- \infty}}^{\infty}\;{\frac{\tan\left( {p_{n}b} \right)}{p_{n}b}\frac{\sin^{2}\left\lbrack {\left( {k_{z} + \frac{2\pi\; n}{L}} \right)\frac{a}{2}} \right\rbrack}{\left\lbrack {\left( {k_{z} + \frac{2\pi\; n}{L}} \right)\frac{a}{2}} \right\rbrack^{2}}}}} +}} \\ {{{\cot\left\lbrack {\frac{\omega}{c}\left( {d - b} \right)} \right\rbrack} = 0},} \end{matrix} & {{Eq}.\mspace{14mu} 29} \end{matrix}$ where p_(n) can be imaginary in which case,

$\begin{matrix} {\frac{\tan\left( {p_{n}b} \right)}{p_{n}b} = {\frac{\tanh\left( {{p_{n}}b} \right)}{{p_{n}}b}.}} & {{Eq}.\mspace{14mu} 30} \end{matrix}$ The mode structures in the double-sided slow-wave structures are qualitatively different from those in the single-sided slow-wave structures, which makes the single-side slow-wave structure suitable for use in a PCRB TWA.

To gain some quantitative understanding of the vacuum dispersion characteristics of electromagnetic wave propagation in the double-side slow-wave structure shown in FIG. 3, consider the first example with the choice of system parameters corresponding to: L=0.24 cm, a/L=0.8, b/L=1.0, d/L=6.0.  Eq. 31 This choice of system parameters can also represent a single-side slow-wave structure that supports only the anti-symmetric modes.

FIGS. 4A-4B show plots of the dispersion functions D_(a)(ω,k) for the anti-symmetric modes as they vary with frequency at several values of phase shifts (or wave numbers). The phase shift is equal to 360°×(k_(z)L/2π). For a given phase shift, the zeros of the dispersion function correspond to the eigenfrequencies of the system. In this example, there is one zero below 20 GHz at least, as shown in FIG. 4A.

FIGS. 5A-5B show plots of the dispersion functions D_(s)(ω,k) for the symmetric modes as they vary with frequency at several values of phase shifts (or wave numbers). The phase shift is equal to 360°×(k_(z)L/2π). For a given phase shift, the zeros of the dispersion function correspond to the eigenfrequencies of the system. In this example, there is one zero below 20 GHz at least, as shown in FIG. 5A.

The eigenfrequencies for the two lowest bands of anti-symmetric modes and the two lowest bands of symmetric modes are plotted as a function of phase shift in FIG. 6.

In the double-sided slow-wave structure, both anti-symmetric and symmetric modes exist. In the phase shift range from 90° to 270°, the anti-symmetric and symmetric modes are nearly degenerate in the first band as well as in the second band, as shown in FIG. 6. In the single-side slow-wave structure, however, only the anti-symmetric modes shown in FIG. 6 exist, and they are well separated in frequency. This sharp contrast in the mode properties makes the single-sided slow-wave structure a better choice for use in the PCRB TWA, because there are no symmetric modes.

The eigenfrequencies of the TWA can also be determined using SUPERFISH code. FIG. 7 shows the comparison between the eigenfrequencies calculated with SUPERFISH and Eq. 27 for the two lowest anti-symmetric modes with the same system parameters as in FIG. 6 or in Eq. 31. Only two lowest modes are plotted in FIG. 7. As can be seen from FIG. 7, the theoretical results and SUPERFISH simulations are in excellent agreement. Therefore, the theory derived for PCRB TWA is applicable for designing an actual device.

The ribbon beam, as shown in FIGS. 1A, is an important component for the efficient operation of the invention. Note the ribbon beam uses a planar wiggler magnetic field.

Let us first consider a ribbon beam in the planar wiggler field B _(w)(x)=−B _(w) [ê _(x) cos h(k _(w) x)cos(k _(w) z)−ê _(z) sin h(k _(w) x)sin(k _(w) z)],  Eq. 32 where B_(w)=constant, k_(w)=2π/λ_(w), and λ_(w) is the wiggler period. Introducing the vector potential

$\begin{matrix} {{A_{w}(x)} = {\frac{m_{e}c^{2}a_{w}}{e}{\cosh\left( {k_{w}x} \right)}{\sin\left( {k_{w}z} \right)}{\hat{e}}_{y}}} & {{Eq}.\mspace{14mu} 33} \end{matrix}$ with the gauge condition B _(w)(x)=∇×A _(w)(x),  Eq. 34 the Hamiltonian for the single-particle motion in cgs units,

$\begin{matrix} {{{H\left( {x,z,P_{x},P_{y},P_{z}} \right)} = {\frac{1}{2m_{e}}\left\lbrack {P_{x}^{2} + \left( {P_{y} + {\frac{e}{c}A_{w}}} \right)^{2} + P_{z}^{2}} \right\rbrack}},} & {{Eq}.\mspace{14mu} 35} \end{matrix}$ can be expanded for P_(y)=mv_(y)−(e/c)A_(w)(x,z)=0 and |k_(w)x|<1 as H(x,z,P _(x) ,P _(y) ,P _(z))≅H ₀(P _(z))+H _(β)(x,P _(x)).  Eq. 36 In Eq. 36, the Hamiltonians

$\begin{matrix} {{{H_{0}\left( P_{z} \right)} = {\frac{P_{z}^{2}}{2m_{e}} + {\frac{1}{4}a_{w}^{2}m_{e}c^{2}}}}{and}} & {{Eq}.\mspace{14mu} 37} \\ {{H_{\beta}\left( {x,P_{x}} \right)} = {\frac{1}{2m_{e}}\left( {P_{x}^{2} + {m_{e}^{2}\omega_{\beta}^{2}x^{2}}} \right)}} & {{Eq}.\mspace{14mu} 38} \end{matrix}$ describe the axial motion and the (transverse) betatron oscillations, respectively, and

$\begin{matrix} {\omega_{\beta} = {{\frac{1}{\sqrt{2}}{ck}_{w}a_{w}} = {\frac{1}{\sqrt{2}}\left( \frac{{eB}_{w}}{m_{e}c} \right)}}} & {{Eq}.\mspace{14mu} 39} \end{matrix}$ is the betatron oscillation frequency.

For P_(y)=0 and |k_(w)x|<1, the equations of motion are

$\begin{matrix} {{\frac{\mathbb{d}z}{\mathbb{d}t} = \frac{P_{z}}{m_{e}}},} & {{Eq}.\mspace{14mu} 40} \\ {{\frac{\mathbb{d}x}{\mathbb{d}t} = \frac{P_{x}}{m_{e}}},} & {{Eq}.\mspace{14mu} 41} \\ {{\frac{\mathbb{d}P_{z}}{\mathbb{d}t} = 0},} & {{Eq}.\mspace{14mu} 42} \\ {\frac{\mathbb{d}P_{x}}{\mathbb{d}t} = {{- m_{e}}\omega_{\beta}^{2}{x.}}} & {{Eq}.\mspace{14mu} 43} \end{matrix}$ Since ν_(z)=β_(b)c=constant,  Eq. 44 the equation of motion for the betatron oscillations can be expressed alternatively as

$\begin{matrix} {{{\frac{\mathbb{d}^{2}x}{\mathbb{d}z^{2}} + {k_{\beta}^{2}x}} = 0},{where}} & {{Eq}.\mspace{14mu} 45} \\ {k_{\beta} = \frac{\omega_{\beta}}{\beta_{b}c}} & {{Eq}.\mspace{14mu} 46} \end{matrix}$ is the betatron wavenumber.

For simplicity, the coupling between the beam envelopes in the x- and y-directions is ignored, and express the beam envelope equation in the x-direction as

$\begin{matrix} {{{\frac{\mathbb{d}^{2}x_{b}}{\mathbb{d}z^{2}} + {k_{\beta}^{2}x_{b}} - \frac{K}{2\left( {x_{b} + y_{b}} \right)} - \frac{ɛ_{x}^{2}}{x_{b}^{3}}} = 0},{where}} & {{Eq}.\mspace{20mu} 47} \\ {{x_{b}^{2} = \left\langle x^{2} \right\rangle},} & {{Eq}.\mspace{14mu} 48} \\ {{y_{b}^{2} = \left\langle y^{2} \right\rangle},} & {{Eq}.\mspace{14mu} 49} \\ {{K = \frac{2I_{b}}{I_{A}\beta_{b}^{2}}},} & {{Eq}.\mspace{14mu} 50} \\ {{\frac{I_{A}}{\beta_{b}} = {\frac{m_{e}c^{3}}{e} \cong {17\mspace{14mu}{kA}}}},} & {{Eq}.\mspace{14mu} 51} \\ {ɛ_{x}^{2} = {{{\left\langle x^{2} \right\rangle\left\langle \left( \frac{\mathbb{d}x}{\mathbb{d}z} \right)^{2} \right\rangle} - \left\langle {x\frac{\mathbb{d}x}{\mathbb{d}z}} \right\rangle^{2}} = {\frac{x_{b}^{2}}{\beta_{b}^{2}c^{2}}{\left( \frac{k_{B}T}{m_{e}} \right).}}}} & {{Eq}.\mspace{14mu} 52} \end{matrix}$ Here, k_(B) is the Boltzmann constant, T is the Kelvin temperature, and <X> denotes

$\begin{matrix} {\left\langle \chi \right\rangle = \frac{\int^{\;}\;{{xf}\ {\mathbb{d}x}{\mathbb{d}y}{\mathbb{d}x^{\prime}}{\mathbb{d}y^{\prime}}}}{\int{f{\mathbb{d}x}{\mathbb{d}y}{\mathbb{d}x^{\prime}}{\mathbb{d}y^{\prime}}}}} & {{Eq}.\mspace{14mu} 53} \end{matrix}$ with x′=dx/dz, y′=dy/dz, and ƒ(x,y,x′,y′) is the electron distribution function.

In the zero-current limit, K=0 and it follows from Eq. 47 that the equilibrium rms beam envelope is given by

$\begin{matrix} {{x_{b} = {\left( \frac{k_{B}T}{m_{e}} \right)^{1/2}\frac{1}{k_{\beta}\beta_{b}c}}},} & {{Eq}.\mspace{14mu} 54} \end{matrix}$ or the wiggler field required for the beam equilibrium amplitude to occur is given by

$\begin{matrix} {{\frac{{eB}_{w}}{m_{e}c} = {\left( \frac{2k_{B}T}{m_{e}} \right)^{1/2}\frac{1}{x_{b}}}},} & {{Eq}.\mspace{14mu} 55} \end{matrix}$ which allows us to calculate B_(w). Note that the wiggler period λ_(w) does not appear in Eq. 55. However, one must demand 2k_(w)x_(b)≦1,  Eq. 56 or λ_(w)≧4πx_(b)  Eq. 57 for the approximations in Eq. 51 to be valid.

For example, taking

$\begin{matrix} {{x_{b} = {0.06\mspace{14mu}{cm}}},{\beta_{b} = 0.08},{{\frac{1}{2}m_{e}\beta_{b}^{2}c^{2}} = {1700\mspace{14mu} e\; V}},{{\frac{1}{2}k_{B}T} = {{\frac{1}{40} \times \frac{1200}{300}} = {0.1\mspace{20mu} e\; V}}},} & {{Eq}.\mspace{14mu} 58} \end{matrix}$ one can obtain from Eq. 55, B_(w)=24.7 G,  Eq. 59 which is easily achievable.

For a finite beam current, K≠0 and it follows from Eq. 47 that the rms beam envelope is given by

$\begin{matrix} {{x_{b} = {\frac{1}{k_{\beta}}\left\lbrack {\frac{K}{2\left( {\xi + 1} \right)} + \frac{k_{B}T}{m_{e}\beta_{b}^{2}c^{2}}} \right\rbrack}^{1/2}},{or}} & {{Eq}.\mspace{14mu} 60} \\ {{\frac{{eB}_{w}}{m_{e}c} = {\frac{\beta_{b}c}{x_{b}}\left\lbrack {\frac{K}{\left( {\xi + 1} \right)} + \frac{2k_{B}T}{m_{e}\beta_{b}^{2}c^{2}}} \right\rbrack}^{1/2}},{where}} & {{Eq}.\mspace{20mu} 61} \\ {\xi = \frac{y_{b}}{x_{b}}} & {{Eq}.\mspace{20mu} 62} \end{matrix}$ is the aspect ratio of the ribbon beam. In this discussion, the value of ξ is fixed.

For example, taking

$\begin{matrix} {{x_{b} = {0.06\mspace{14mu}{cm}}},{\beta_{b} = 0.08},{{\frac{1}{2}m_{e}\beta_{b}^{2}c^{2}} = {1700\mspace{14mu} e\; V}},{{\frac{1}{2}k_{B}T} = {{\frac{1}{40} \times \frac{1200}{300}} = {0.1\mspace{14mu} e\; V}}},{I_{b} = {0.076\mspace{14mu} A}},{\xi = 10},} & {{Eq}.\mspace{20mu} 63} \end{matrix}$ one obtains from Eq. 63, B_(w)=93.9 G,  Eq.64 which is still easily achievable.

Now, a ribbon electron beam 46 interacts with a single-sided PC slow-wave structure 48, similar to structure 38 of FIG. 2, shown in FIG. 8. The equilibrium of the beam 46 is described by V=V_(b)ê_(z)=β_(b)cê_(z),  Eq. 65

$\begin{matrix} {{J = {{{- e}\;\sigma_{b}V_{b}{\delta\left( {x - h} \right)}{\hat{e}}_{z}} = {{- \frac{\mathbb{d}I_{b}}{\mathbb{d}y}}{\delta\left( {x - h} \right)}{\hat{e}}_{z}}}},} & {{Eq}.\mspace{14mu} 66} \end{matrix}$ where V_(b)=constant is the equilibrium beam velocity, σ_(b) is the surface number density of the electrons in the ribbon beam, the sheet x=h specifies the transverse displacement of the beam, and

$\begin{matrix} {I_{b} = {\int_{y = {{- 2}y_{b}}}^{2y_{b}}{\int_{x = 0}^{b}{\frac{\mathbb{d}I_{b}}{\mathbb{d}y}{\delta\left( {x - h} \right)}\ {\mathbb{d}x}\ {\mathbb{d}y}}}}} & {{Eq}.\mspace{14mu} 67} \end{matrix}$ is the beam current, and y_(b) is the rms width of the ribbon beam.

In this analysis, the variations are ignored in the y-direction, and treat the stability of the ribbon beam 46 as a two-dimensional problem. The linearized cold-fluid equations are

$\begin{matrix} {{{\left( {\frac{\partial\;}{\partial t} + {V_{b}\frac{\partial\;}{\partial z}}} \right)\delta\; V_{z}} = {{{- \frac{e}{m}}\delta\; E_{z}}❘_{x = h}}},} & {{Eq}.\mspace{20mu} 68} \\ {{{{\left( {\frac{\partial}{\partial t} + {V_{b}\frac{\partial\;}{\partial z}}} \right){\delta\sigma}} + {\sigma_{b}\frac{\partial\;}{\partial z}\delta\; V_{z}}} = 0},} & {{Eq}.\mspace{14mu} 69} \\ \begin{matrix} {{\left( {\frac{\partial^{2}}{\partial x^{2}} + \frac{\partial^{2}}{\partial z^{2}} - {\frac{1}{c^{2}}\frac{\partial^{2}}{\partial t^{2}}}} \right)\delta\; E_{z}} = {{\frac{4\pi}{c^{2}}\frac{{\partial\delta}\; J_{z}}{\partial t}} + {4\pi\frac{\partial{\delta\rho}}{\partial z}}}} \\ {= {{- 4}\pi\; e\left\{ {\frac{1}{c^{2}}\frac{\partial\;}{\partial t}\left( {{V_{b}{\delta\sigma}} +} \right.} \right.}} \\ {{\left. {\left. {\sigma_{b}\delta\; V_{z}} \right) + {\frac{\partial\;}{\partial z}{\delta\sigma}}} \right\}{\delta\left( {x - h} \right)}},} \end{matrix} & {{Eq}.\mspace{14mu} 70} \end{matrix}$ where the charge and current density perturbations δρ(x,z,t) and δJ_(z)(x,z,t) are defined as δρ(x,z,t)=−eδσ(z,t)δ(x−h),  Eq. 71 δJ _(z)(x,z,t)=−eδσ(z,t)δ(x−h)V _(b) −eσ _(b)δ(x−h)δV _(z)(z,t).  Eq. 72 Expressing all perturbations as

$\begin{matrix} {{{{\delta\psi}\left( {x,z,t} \right)} = {{\mathbb{e}}^{{\mathbb{i}}{({{k_{z}z} - {\omega\; t}})}}{\sum\limits_{n = {- \infty}}^{\infty}\;{{{\delta\psi}_{n}(x)}{\mathbb{e}}^{{\mathbb{i}2\pi}\;{{nz}/L}}}}}},{or}} & {{Eq}.\mspace{14mu} 73} \\ {{{\delta\;{V_{z}\left( {z,t} \right)}} = {{\mathbb{e}}^{{\mathbb{i}}{({{k_{z}z} - {\omega\; t}})}}{\sum\limits_{n = {- \infty}}^{\infty}\;{\delta\; V_{2n}{\mathbb{e}}^{{\mathbb{i}2\pi}\;{{nz}/L}}}}}},} & {{Eq}.\mspace{14mu} 74} \\ {{{{\delta\sigma}\left( {z,t} \right)} = {{\mathbb{e}}^{{\mathbb{i}}{({k,{z - {\omega\; t}}})}}{\sum\limits_{n = {- \infty}}^{\infty}\;{{\delta\sigma}_{n}{\mathbb{e}}^{{\mathbb{i}2\pi}\;{{nz}/L}}}}}},} & {{Eq}.\mspace{14mu} 75} \\ {{{\delta\;{E_{z}\left( {x,z,t} \right)}} = {{\mathbb{e}}^{{\mathbb{i}}{({k,{z - {\omega\; t}}})}}{\sum\limits_{n = {- \infty}}^{\infty}\;{\delta\;{E_{2n}(x)}{\mathbb{e}}^{{\mathbb{i}2\pi}\;{{nz}/L}}}}}},} & {{Eq}.\mspace{14mu} 76} \end{matrix}$ one can rewrite the linearized cold-fluid equations 69-70 as

$\begin{matrix} {{{{i\left( {\omega - {k_{n}V_{b}}} \right)}\mspace{11mu}\delta\; V_{zn}} = {\frac{e}{m}\delta\;{E_{sn}(h)}}},} & {{Eq}.\mspace{14mu} 77} \end{matrix}$ (ω−k _(n) V _(b))δσ_(n) −k _(n)σ_(b) δV _(zn)=0,  Eq. 78

$\begin{matrix} {{{\left( {\frac{\mathbb{d}^{2}}{\mathbb{d}x^{2}} + p_{n}^{2}} \right)\delta\;{E_{zn}(x)}} = {4\pi\;{{ie}\left( {{\frac{\omega\; V_{b}}{c^{2}}{\delta\sigma}_{n}} + {\frac{{\omega\sigma}_{b}}{c^{2}}\delta\; V_{zn}} - {k_{n}{\delta\sigma}_{n}}} \right)}{\delta\left( {x - h} \right)}}},} & {{Eq}.\mspace{14mu} 79} \end{matrix}$ where k_(n) and p_(n) ² are defined as

$\begin{matrix} {{k_{n} = {k_{z} + \frac{2\pi\; n}{L}}},} & {{Eq}.\mspace{14mu} 80} \\ {p_{n}^{2} = {\frac{\omega^{2}}{c^{2}} - {k_{n}^{2}.}}} & {{Eq}.\mspace{14mu} 81} \end{matrix}$ Solving for δV_(zn) and δσ_(n) from Eqs. 67 and 68 in terms of δE_(zn)(h) yields

$\begin{matrix} {{{\delta\; V_{zn}} = {{- \frac{ie}{m}}\frac{\delta\;{E_{zn}(h)}}{\left( {\omega - {k_{n}V_{b}}} \right)}}},} & {{Eq}.\mspace{14mu} 82} \\ {{\delta\sigma}_{n} = {{- \frac{ie}{m}}{\frac{k_{n}\sigma_{b}\delta\;{E_{2n}(h)}}{\left( {\omega - {k_{n}V_{b}}} \right)^{2}}.}}} & {{Eq}.\mspace{14mu} 83} \end{matrix}$ Substituting Eqs. 82 and 83 into Eq. 79, one obtains

$\begin{matrix} {{\left( {\frac{\mathbb{d}^{2}}{\mathbb{d}x^{2}} + p_{n}^{2}} \right)\delta\;{E_{zn}(x)}} = {\frac{4{\pi\mathbb{e}}^{2}\sigma_{b}}{m}\frac{p_{n}^{2}\delta\;{E_{zn}(h)}}{\left( {\omega - {k_{n}V_{b}}} \right)^{2}}{{\delta\left( {x - h} \right)}.}}} & {{Eq}.\mspace{14mu} 84} \end{matrix}$ The solutions to Eq. 84 are

$\begin{matrix} {{\delta\;{E_{zn}(x)}} = \left\{ \begin{matrix} {{{A_{n}^{<}\mspace{14mu}{\sin\left( {p_{n}x} \right)}} + {B_{n}^{<}{\cos\left( {p_{n}x} \right)}}},} & {{0 < x < h},} \\ {{{A_{n}^{>}{\sin\left( {p_{n}x} \right)}} + {B_{n}^{>}{\cos\left( {p_{n}x} \right)}}},} & {{h < x < b},} \end{matrix} \right.} & {{Eq}.\mspace{20mu} 85} \end{matrix}$ where A_(n) ^(<), B_(n) ^(<), A_(n) ^(>) and B_(n) ^(>) are constants to be determined by the boundary conditions at x=0, x=h, and x=b.

The coefficients A_(n) ^(<), B_(n) ^(<), A_(n) ^(>) and B_(n) ^(>) are determined in the following three steps. First, note the electric field δE_(z)(x,z,t) vanishes at x=0. This requires δE _(zn)(0)=0.  Eq. 86 Therefore, one must have B _(n) ^(<)=0.  Eq. 87

Second, the coefficients A_(n) ^(>) and B_(n) ^(>) are expressed in terms of A_(n) ^(<), using both the continuity of the axial electric field at x=h, for example, A _(n) ^(>)sin(p _(n) h)+B _(n) ^(>)cos(p _(n) h)−A _(n) ^(<)sin(p _(n) h)=0,  Eq. 88 and the relation

$\begin{matrix} {{{\frac{{\mathbb{d}\delta}\; E_{zn}}{\mathbb{d}x}}_{x = {h - 0}}^{h + 0} = {\frac{4{\pi\mathbb{e}}^{2}\sigma_{b}}{m}\frac{p_{n}^{2}}{\left( {\omega - {k_{n}V_{b}}} \right)^{2}}A_{n}^{<}{\sin\left( {p_{n}h} \right)}}}{or}} & {{Eq}.\mspace{14mu} 89} \\ {{{A_{n}^{>}{\cos\left( {p_{n}h} \right)}} - {B_{n}^{>}{\sin\left( {p_{n}h} \right)}} - {A_{n}^{<}{\cos\left( {p_{n}h} \right)}}} = {\frac{4{\pi\mathbb{e}}^{2}\sigma_{b}}{m}\frac{p_{n}}{\left( {\omega - {k_{n}V_{b}}} \right)^{2}}A_{n}^{<}\mspace{11mu}{{\sin\left( {p_{n}h} \right)}.}}} & {{Eq}.\mspace{14mu} 90} \end{matrix}$ Solving Eqs. 88 and 90 for A_(n) ^(>) and B_(n) ^(>) gives

$\begin{matrix} {{\frac{A_{n}^{>}}{A_{n}^{<}} = {1 + {{\alpha_{n}\left( {\omega,k} \right)}{\cos\left( {p_{n}h} \right)}}}},} & {{Eq}.\mspace{14mu} 91} \\ {{\frac{B_{n}^{>}}{A_{n}^{<}} = {{- {\alpha_{n}\left( {\omega,k} \right)}}{\sin\left( {p_{n}h} \right)}}},} & {{Eq}.\mspace{14mu} 92} \end{matrix}$ where the function α_(n)(ω,k) is defined by

$\begin{matrix} {{\alpha_{n}\left( {\omega,k} \right)} = {\frac{4\pi\;{\mathbb{e}}^{2}\sigma_{b}}{m_{e}}{\frac{p_{n}\mspace{14mu}{\sin\left( {p_{n}h} \right)}}{\left( {\omega - {k_{n}V_{b}}} \right)^{2}}.}}} & {{Eq}.\mspace{14mu} 93} \end{matrix}$ Furthermore, the electric and magnetic fields are expressed in the region h<x<b as

$\begin{matrix} {{{\delta\;{E_{z}\left( {x,z,t} \right)}} = {\sum\limits_{n = {- \infty}}^{\infty}\;{\left\lbrack {{A_{n}^{>}{\sin\left( {p_{n}x} \right)}} + {B_{n}^{>}{\cos\left( {p_{n}x} \right)}}} \right\rbrack{\mathbb{e}}^{{\mathbb{i}}{({{k_{n}z} - {\omega\; t}})}}}}},} & {{Eq}.\mspace{14mu} 94} \\ {{{\delta\;{E_{x}\left( {x,z,t} \right)}} = {\frac{c}{i\;\omega}\frac{\partial\;}{\partial z}\delta\;{B_{y}\left( {x,z,t} \right)}}},} & {{Eq}.\mspace{14mu} 95} \\ {{\delta\;{B_{y}\left( {x,z,t} \right)}} = {\frac{{\mathbb{i}}\;\omega}{{cp}_{n}}{\sum\limits_{n = {- \infty}}^{\infty}{\left\lbrack {{A_{n}^{>}{\cos\left( {p_{n}x} \right)}} - {B_{n}^{>}{\cos\left( {p_{n}x} \right)}}} \right\rbrack{{\mathbb{e}}^{{\mathbb{i}}{({{k_{n}z} - {\omega\; t}})}}.}}}}} & {{Eq}.\mspace{14mu} 96} \end{matrix}$

Third, a relationship is derived between the coefficient A_(n) ^(<) and the electric field amplitude E₀ in the corrugated-vane region, using the approximations in Eqs. 11-13,

$\begin{matrix} {{{B_{y}\left( {x,z,t} \right)} \cong {{iE}_{0}{\mathbb{e}}^{{{\mathbb{i}}{({{k_{z}{sL}} - {\omega\; t}})}}\;}{\cos\left\lbrack {\frac{\omega}{c}\left( {x - d} \right)} \right\rbrack}}},} & {{Eq}.\mspace{14mu} 97} \end{matrix}$ E _(x)(x,z,t)≅0,  Eq. 98

$\begin{matrix} {{E_{z}\left( {x,z,t} \right)} \cong {E_{0}{\mathbb{e}}^{{\mathbb{i}}{({{k_{z}{sL}} - {\omega t}})}}{\sin\left\lbrack {\frac{\omega}{c}\left( {x - d} \right)} \right\rbrack}}} & {{Eq}.\mspace{14mu} 99} \end{matrix}$ for |z−sL|<a/2 and b<x<d. Here, s=0, ±1, ±2, . . . This gives

$\begin{matrix} {{{\left\lbrack {{\frac{A_{n}^{>}}{A_{n}^{<}}{\sin\left( {p_{n}b} \right)}} + {\frac{B_{n}^{>}}{A_{n}^{<}}{\cos\left( {p_{n}b} \right)}}} \right\rbrack A_{n}^{<}} = {\frac{{aE}_{0}}{L}{\sin\left\lbrack {\frac{\omega}{c}\left( {b - d} \right)} \right\rbrack}\frac{\sin\left\lbrack {\left( {k_{z} + \frac{2\pi\; n}{L}} \right)\frac{a}{2}} \right\rbrack}{\left( {k_{z} + \frac{2\pi\; n}{L}} \right)\frac{a}{2}}}},} & {{Eq}.\mspace{14mu} 100} \end{matrix}$ which relates A_(n) ^(<) and E₀.

Using the expressions for the electric and magnetic fields in Eqs. 94-99, the continuity of the axial electric field at x=b, and Eq. 100, it is readily shown that the loaded admittances A_(L) ^(±)(ω,k_(z)), defined in the same manner as the unloaded admittances A^(±)(ω,k_(z)) in Eq. 23, are given by

$\begin{matrix} {{{A_{L}^{+}\left( {\omega,k_{z}} \right)} = {{A^{+}\left( {\omega,k_{z}} \right)} = {\frac{1}{ia}{{\cot\left\lbrack {\frac{\omega}{c}\left( {d - b} \right)} \right\rbrack} \cdot \frac{\left( \frac{k_{z}a}{2} \right)^{2}}{\sin^{2}\left( \frac{k_{z}a}{2} \right)}}}}}{and}} & {{Eq}.\mspace{14mu} 101} \\ {{{{A_{L}^{-}\left( {\omega,k_{z}} \right)} = {{\left( \frac{{\mathbb{i}}\;\omega\; b}{cL} \right){\frac{\frac{k_{z}a}{2}}{\sin\left( \frac{k_{z}a}{2} \right)}}^{2} \times {\sum\limits_{n = {- \infty}}^{\infty}\;{\left( \frac{1}{p_{n}b} \right)\frac{{\frac{A_{n}^{>}}{A_{n}^{<}}{\cos\left( {p_{n}b} \right)}} - {\frac{B_{n}^{>}}{A_{n}^{<}}{\sin\left( {p_{n}b} \right)}}}{{\frac{A_{n}^{>}}{A_{n}^{<}}{\sin\left( {p_{n}b} \right)}} + {\frac{B_{n}^{>}}{A_{n}^{<}}{\cos\left( {p_{n}b} \right)}}}}}}❘\frac{\sin\left\lbrack {\left( {k_{z} + \frac{2\pi\; n}{L}} \right)\frac{a}{2}} \right\rbrack}{\left( {k_{z} + \frac{2\pi\; n}{L}} \right)\frac{a}{2}}}}}^{2}.} & {{Eq}.\mspace{14mu} 102} \end{matrix}$ Setting A _(L) ⁺(ω,k _(z))=A _(L) ⁻(ω,k _(z)),  Eq. 103 one obtains the loaded dispersion relation

$\begin{matrix} {{{D_{a}^{L}\left( {\omega,k_{z}} \right)} \equiv {{\cot\left\lbrack {\frac{\omega}{c}\left( {d - b} \right)} \right\rbrack} + {\frac{\omega\; a\; b}{cL}{\overset{\infty}{\sum\limits_{n = {- \infty}}}{\left( \frac{1}{p_{n}b} \right)\frac{{\cos\left( {p_{n}b} \right)} + {\alpha_{n}{\cos\left\lbrack {p_{n}\left( {b - h} \right)} \right\rbrack}}}{{\sin\left( {p_{n}b} \right)} + {\alpha_{n}{\sin\left\lbrack {p_{n}\left( {b - h} \right)} \right\rbrack}}}{\frac{\sin\left( {k_{n}\frac{a}{2}} \right)}{k_{n}\frac{a}{2}}}^{2}}}}}} = 0} & {{Eq}.\mspace{14mu} 104} \end{matrix}$ where use has been made of Eqs. 91-93, and the functions k_(n)=k_(n)(k_(z)), p_(n)=p_(n)(ω,k_(z)) and α_(n)=α_(n)(ω,k_(z)) are defined in Eqs. 80, 81 and 93, respectively.

When I_(b)=0 (or α_(n)=0), the dispersion relation (Eq. 104) reduces to the vacuum dispersion relation (Eq. 91) as expected.

In the Compton regime, |α_(n)|<<1 and one may make the approximation

                                        Eq.  105 ${\frac{{\cos\left( {p_{n}b} \right)} + {\alpha_{n}{\cos\left\lbrack {p_{n}\left( {b - h} \right)} \right\rbrack}}}{{\sin\left( {p_{n}b} \right)} + {\alpha_{n}{\sin\left\lbrack {p_{n}\left( {b - h} \right)} \right\rbrack}}} \cong {{\cot\left( {p_{n}b} \right)} + {\alpha_{n}\frac{\sin\left( {p_{n}h} \right)}{\sin^{2}\left( {p_{n}b} \right)}}}}\mspace{304mu} = {{\cot\left( {p_{n}b} \right)} + {\frac{4\pi\; e^{2}\sigma_{b}}{m_{e}}\frac{p_{n}}{\left( {\omega - {k_{n}V_{b}}} \right)^{2}}\frac{\sin^{2}\left( {p_{n}h} \right)}{\sin^{2}\left( {p_{n}b} \right)}}}$ and express the loaded dispersion relation approximately as

$\begin{matrix} {{{D_{a}\left( {\omega,k_{z}} \right)} = {\frac{\omega\; a}{c\; L}\left( \frac{4\pi\; e^{2}\sigma_{b}}{m_{e}} \right){\overset{\infty}{\sum\limits_{n = {- \infty}}}{\frac{1}{\left( {\omega - {k_{n}V_{b}}} \right)^{2}}\frac{\sin^{2}\left( {p_{n}h} \right)}{\sin^{2}\left( {p_{n}b} \right)}{\frac{\sin\left( {k_{n}\frac{a}{2}} \right)}{k_{n}\frac{a}{2}}}^{2}}}}},} & {{Eq}.\mspace{14mu} 106} \end{matrix}$ where D_(a)(ω,k_(z)) is the vacuum dispersion function defined in Eq. (3.1.27). Furthermore, if only one term on the right-hand side of Eq. (5.2.2) dominates, say n=m, then one can further approximate the loaded dispersion relation as D _(a)(ω,k _(z))(ω−k _(m) V _(b))²={tilde over (ε)}_(m)(ω,k _(z)),  Eq. 107 where the coupling parameter {tilde over (ε)}_(m)(ω,k_(z)) is defined by

$\begin{matrix} {{{\overset{\sim}{ɛ}}_{m}\left( {\omega,k_{z}} \right)} = {\frac{\omega\; a}{c\; L}\left( \frac{4\pi\; e^{2}\sigma_{b}}{m_{e}} \right)\frac{\sin^{2}\left( {p_{m}h} \right)}{\sin^{2}\left( {p_{m}b} \right)}{{\frac{\sin\left( {k_{m}\frac{a}{2}} \right\rbrack}{k_{m}\frac{a}{2}}}^{2}.}}} & {{Eq}.\mspace{14mu} 108} \end{matrix}$

To estimate the linear gain and bandwidth in the Compton regime, let (ω_(c),k_(c)) denote an intersection point of D _(a)(ω,k _(z))=0  Eq. 109 and ω−(k _(z)+2πm/L)V _(b)=0  Eq. 110 in the ω versus k_(z) diagram. Expanding D_(a)(ω,k_(z)) about the point (ω,k_(z))=(ω_(c),k_(c)), i.e.,

$\begin{matrix} {{{D_{a}\left( {\omega,k_{z}} \right)} \cong {{\left( \frac{\partial D_{a}}{\partial\omega} \right)_{({\omega_{c},k_{c}})}\left( {\omega - \omega_{c}} \right)} + {\left( \frac{\partial D_{a}}{\partial k_{z}} \right)_{({\omega_{c},k_{c}})}\left( {k_{z} - k_{c}} \right)}}}\mspace{104mu} = {\left( \frac{\partial D_{a}}{\partial\omega} \right)_{({\omega_{c},k_{c}})}\left\lbrack {\omega - \omega_{c} - {\upsilon_{g}\left( {k_{z} - k_{c}} \right)}} \right\rbrack}} & {{Eq}.\mspace{14mu} 111} \end{matrix}$ with ν_(g) being the group velocity, and introducing the rescaled coupling parameter

$\begin{matrix} {{ɛ_{m} = {\left( \frac{\partial D_{a}}{\partial\omega} \right)_{({\omega_{c},k_{c}})}^{- 1}{{\overset{\sim}{ɛ}}_{m}\left( {\omega_{c},k_{c}} \right)}}},} & {{Eq}.\mspace{14mu} 112} \end{matrix}$ one can express the loaded dispersion relation in the following simplified form [ω−ω_(c)−ν_(g)(k _(z) −k _(c))][ω−k _(z) V _(b)−(2πm/L)V _(b)]²=ε_(m).  Eq. 113 Because Eq. 113 is cubic in either ω or k_(z), it can be solved analytically.

The frequency shift δω and the detuning parameter ΔΩ_(m) is defined as δω=ω−k _(z) V _(b)−(2πm/L)V _(b),  Eq. 114 ΔΩ_(m)=ω_(c)+ν_(g)(k _(z) −k _(c))−k _(z) V _(b)−(2πm/L)V _(b).  Eq. 115 The loaded dispersion relation becomes (δω)²(δω−ΔΩ_(m))=ε_(m).  Eq. 116 Since {tilde over (ε)}_(m)>0 and (∂D_(a)/∂ω)<0 (see FIG. 3), where ε_(m)<0. Consequently, Eq. 116 yields the maximum temporal growth rate

$\begin{matrix} {{{{Im}\;{\delta\omega}}}_{\max} = {\frac{\sqrt{3}}{2}{ɛ_{m}}^{\frac{1}{3}}}} & {{Eq}.\mspace{14mu} 117} \end{matrix}$ at ΔΩ_(m)=0. Since ε_(m)∝I_(b), the scaling relation |Imδω|_(max)∝I_(b) ^(1/3)  Eq. 118 holds in the Compton regime.

In the Raman regime, one must treat the space-charge term in Eq. 104 carefully. In this case, the Eq. 104 is expressed as

$\begin{matrix} {{D_{a}\left( {\omega,k_{z}} \right)} = {\frac{\omega\; a\; b}{cL}{\overset{\infty}{\sum\limits_{n = {- \infty}}}{\left( \frac{1}{p_{n}b} \right)\left\{ {\frac{{\cos\left( {p_{n}b} \right)} + {\alpha_{n}{\cos\left\lbrack {p_{n}\left( {b - h} \right)} \right\rbrack}}}{{\sin\left( {p_{n}b} \right)} + {\alpha_{n}{\sin\left\lbrack {p_{n}\left( {b - h} \right)} \right\rbrack}}} - {\cot\left( {p_{n}b} \right)}} \right\}{\frac{\sin\left( {k_{n}\frac{a}{2}} \right\rbrack}{k_{n}\frac{a}{2}}}^{2}}}}} & {{Eq}.\mspace{14mu} 119} \end{matrix}$ where D_(a)(ω,k_(z)) is the vacuum dispersion function defined in Eq. 27. Substituting Eq. 93 into Eq. 119, and assuming that only one term on the right-hand side of Eq. 106 dominates, say n=m, D _(a)(ω,k _(z))└(ω−k _(m) V _(b))²−(QC)_(m) ²(ω,k _(z))┘={tilde over (ε)}_(m)(ω,k _(z)),  Eq. 120 where the coupling parameter {tilde over (ε)}_(m)(ω,k_(z)) is defined in Eq. 108, and the space-charge parameter (QC)_(m) ² is defined by

$\begin{matrix} {{({QC})_{m}^{2}\left( {\omega,k_{z}} \right)} = {{- \left( \frac{4\pi\; e^{2}\sigma_{b}}{m_{e}} \right)}{\frac{p_{m}{\sin\left( {p_{m}h} \right)}{\sin\left\lbrack {p_{m}\left( {b - h} \right)} \right\rbrack}}{\sin^{2}\left( {p_{m}b} \right)}.}}} & {{Eq}.\mspace{14mu} 121} \end{matrix}$ Typically, the space-charge parameter (QC)_(m) ²(ω,k_(z)) is positive in the regime of interest.

To estimate the linear gain and bandwidth in the Raman regime, let (ω_(c),k_(c)) denote an intersection point of D _(a)(ω,k _(z))=0  Eq. 122 and ω−(k _(z)+2πm/L)V _(b)−(QC)_(m)(ω,k _(z))=0  Eq. 123 in the ω versus k_(z) diagram. Making use of the expansion in Eq. 111, one can express the loaded dispersion relation in the following simplified form

$\begin{matrix} {{{\left\lbrack {\omega - \omega_{c} - {\upsilon_{g}\left( {k_{2} - k_{c}} \right)}} \right\rbrack\left\lbrack {\omega - {k_{z}V_{b}} - {\left( {2\pi\;{m/L}} \right)V_{b}} - {({QC})_{m}\left( {\omega_{c},k_{c}} \right)}} \right\rbrack} = \frac{ɛ_{m}}{2({QC})_{m}\left( {\omega_{c},k_{c}} \right)}},} & {{Eq}.\mspace{14mu} 124} \end{matrix}$ where ε_(m) is defined in Eq. 112.

Following an earlier analysis, the frequency shift δω and the detuning parameter ΔΩ_(m) is defined as δω=ω−k _(z) V _(b)−(2πm/L)V _(b)−(QC)_(m)(ω_(c) ,k _(c)),  Eq. 125 ΔΩ_(m)=ω_(c)+ν_(g)(k _(z) −k _(c))−k _(z) V _(b)−(2πm/L)V _(b)−(QC)_(m)(ω_(c) ,k _(c).  Eq. 126 The loaded dispersion relation becomes

$\begin{matrix} {{{\delta\omega}\left( {{\delta\omega} - {\Delta\Omega}_{m}} \right)} = {\frac{ɛ_{m}}{2({QC})_{m}\left( {\omega_{c},k_{c}} \right)}.}} & {{Eq}.\mspace{14mu} 127} \end{matrix}$ If ε_(m)/(QC)_(m)(ω_(c),k_(c))<0 the system is unstable, and the maximum temporal growth rate is given by

$\begin{matrix} {{{{Im}\mspace{11mu}{\delta\omega}}}_{\max} = {\frac{ɛ_{m}}{2({QC})_{m}\left( {\omega_{c},k_{c}} \right)}}^{1/2}} & {{Eq}.\mspace{14mu} 128} \end{matrix}$ at ΔΩ_(m)=0. Since ε_(m)/(QC)_(m)(ω_(c),k_(c))∝I_(b) ^(1/2), the scaling relation |Imδω|_(max)∝I_(b) ^(1/4)  Eq. 129 holds in the Raman regime.

The dispersion relation in Eq. 104 can be solved numerically using Newton's method to obtain the linear gain. For a real value of the wavenumber k_(z), the temporal linear growth rate ω_(i)=Im(ω)>0 can be obtained from the complex ω that solves Eq. 104. On the other hand, for a real value of the angular frequency ω, the spatial linear growth rate k_(i)=−Im(k_(z))>0 can be obtained from the complex wavenumber k_(z) that solves Eq. 104.

FIGS. 9A-9B show plots of the temporal linear growth rate ω_(i)=Im(ω) of the lowest anti-symmetric mode as a function of the wavenumber k_(z) for a/L=0.8, b/L=1.0, d/L=6.0, h/L=0.8, β_(b)=0.08 and two cases (a) ω_(pb) ²L²/c²=7.0×10⁻⁵ and (b) ω_(pb) ²L²/c²=7.0×10⁻⁷, where ω_(pb) ²L²/c²=4πe²σ_(b)L/c²m_(e)=πI_(b)L/β_(b)I_(c)y_(b) with I_(c)=17 kA. The parameters for the RF circuit in FIGS. 9A-9B are identical to those used in FIGS. 6 and 7. As seen from FIGS. 9A and 9B, each curve has a peak near ω≈k_(z)V_(b), which represents a stronger beam-wave interaction. For the purpose of illustration, the beam parameters in these examples correspond to backward-wave oscillators not traveling-wave amplifiers. For this RF circuit, a more energetic electron beam is required in order to make an amplifier. The high-current case shown in FIG. 9A indicates a temporal growth rate in the Raman regime (i.e., with an upright bell-shaped curve), whereas the low-current case shown in FIG. 9B implies a temporal growth rate in the Compton regime (i.e., with a tilted bell-shaped curve).

FIGS. 10A-10B show plots of the spatial linear growth rate k_(i)=−Im(k_(z)) of the lowest anti-symmetric mode as a function of angular frequency ω for a/L=0.8, b/L=1.0, d/L=6.0, h/L=0.8, β_(b)=0.08. FIG. 9A shows the case ω_(pb) ²L²/c²=7.0×10⁻⁵ and FIG. 9B shows the case ω_(pb) ²L²/c²=7.0×10⁻⁷. These parameters are identical to those in FIGS. 8A-8B.

FIG. 11 shows the maximum temporal growth ω_(i,max)=[Im(ω)]_(max) of the lowest anti-symmetric mode as a function of normalized beam current ω_(pb) ²L²/c²=4πe²σ_(b)L/c²m_(e)=πI_(b)L/β_(b)I_(c)y _(b) for a/L=0.8, b/L=1.0, d/L=6.0, h/L=0.8, and β_(b)=0.08.

FIG. 12 shows the maximum spatial growth k_(i,max)=−[Im(k_(z))]_(max) of the lowest anti-symmetric mode as a function of normalized beam current ω_(pb) ²L²/c²=4πe²σ_(b)L/c²m_(e)=πI_(b)L/β_(b)I_(c)y_(b) for the same parameters as in FIGS. 10A-10B. FIGS. 10A-10B and 11 show that both peak temporals and spatial growth rates increase with the beam current.

It is important to suppress any potential unwanted modes in a microwave amplifier. This is true for the PCRB TWA. In the PCRB TWA, two techniques are used to suppress unwanted modes.

One technique is use of a single-sided slow-wave structure instead of a double-sided slow-wave structure, which eliminates the symmetric modes in the operating band and higher frequency bands.

The other technique is use of photonic crystals. Typically, photonic crystals include periodic metallic structures (e.g., periodic metal rods) or periodic dielectric (e.g., periodic dielectric layers, rods or spheres) or a combination of periodic metallic and dielectric structures. They can be one-, two-, or three-dimensional.

When designed properly, a photonic crystal acts as a frequency-selective and/or mode-selective filter, which keeps the desired operating mode in the amplifier, and at same time, allows other modes, especially unwanted modes, to escape from the amplifier. In other words, the photonic crystal effectively damps the unwanted modes. The effectiveness of photonic crystals in both frequency selection and mode selection were demonstrated in an oscillator operating at high-frequencies and using an oversized cavity with its characteristic size greater than the wavelength, but it remains to be seen in amplifier configurations, especially for transverse size less than the wavelength.

As an example, the dispersion characteristics of wave propagation in photonic crystals can be calculated using the latest Photonic Band Gap Structure Simulator (PBGSS) code developed at MIT. Shown in FIGS. 13A-13B is an example of the dispersion characteristics in a two-dimensional dielectric photonic crystal, where FIG. 13A shows the TM modes and FIG. 13B shows the TE modes, with the rod radius r₀ and the lattice period p, as calculated using the real-space finite-difference (RSFD) and Fourier transform (FT) methods. In this example, there is only a partial band gap at the X point. Here, r₀/p=0.2, ε₁=3.0 for the rods, ε₂=8.9 for the background, and a mesh of 21×21 cells and 13×13 plane waves are used in the PBGSS RSFD and FT calculations, respectively.

The detailed concept design of the PCRB TWA for 3G wireless base stations. will focus on the frequency range from 1920 to 1980 MHz, which is used the initial rollout of 3G wireless network. The PCRB TWA is a 200W, 1950 MHz, 3% bandwidth structure. The parameters and design results are summarized in Table 1 and a cross-sectional view of the amplifier beam tunnel is shown in FIG. 14.

TABLE 1 Circuit Parameters: Structure Single-sided slow-wave structure Axial Period L = 0.478 cm Slot Width a = 0.382 cm Slot Depth d = 4.20 cm Tunnel Height b = 0.478 cm Tunnel Width w = 1.6 cm Beam Parameters: Current I_(b) = 0.11 A Energy ${\frac{1}{2}m_{e}\beta_{b}^{2}c^{2}} = {2290\mspace{14mu}{eV}\mspace{14mu}\left( {\beta_{b} = 0.0946} \right)}$ Displacement h = 0.378 cm Full Height ${4x_{b}} = {0.12\mspace{14mu}{cm}\mspace{20mu}\left( {\frac{2x_{b}}{d - h} = {63\%}} \right)}$ Full Width 4_(Yb) = 1.2 cm Beam Temperature ${\frac{1}{2}k_{B}T} = {0.1\mspace{14mu}{eV}}$ $\begin{matrix} \left\lbrack {{Calculated}\mspace{14mu}{Coupling}\mspace{14mu}{Constant}} \right. \\ {{\frac{\omega_{{pb},{eff}}^{2}L^{2}}{c^{2}} \equiv \frac{4{\pi e}^{2}\sigma_{b}L}{m_{e}c^{2}}} = {\frac{{\pi I}_{b}L}{I_{A}y_{b}} = {3.4 \times 10^{- 4}\mspace{25mu}{with}}}} \\ \left. {{\frac{I_{A}}{\beta_{b}} \equiv \frac{m_{e}c^{3}}{e}} = {{17\mspace{14mu}{kA}\mspace{14mu}{and}\mspace{14mu} I_{b}} = {4e\;\sigma_{b}y_{b}\beta_{b}c}}} \right\rbrack \end{matrix}\quad$ Wiggler Field Parameters: B_(w) = 207 G λ_(w) = 2L = 0.956 cm Results of Gain Calculation: Center Frequency f = 1950 MHz Full Bandwidth Δf/f = 3.1% Gain G = 1.6−3.1 dB/period(3.4-6.6 dB/cm) Efficiency η = 80% Photonic Crystal: Structure Dielectric square lattice Band Gap Width 100 MHz

Shown in FIG. 15 is a plot of the eigenfrequency versus the phase shift for the two lowest bands in the 200 W, 1950 MHz, 3% bandwidth PCRB TWA, as calculated analytically using Eq. 27 with the RF circuit parameters listed in Table 1.

The ribbon electron beam is designed to interact with the lowest band at about 120° phase shift to achieve RF signal amplification. Using the loaded dispersion relation in Eq. 104 and the parameters listed in Table 1, the complex wavenumber k_(z) is calculated using the GAIN code. The results are summarized in FIGS. 16-19.

Shown in FIGS. 16 and 17 are plots of the normalized real and imaginary wavenumbers Rek_(z)L and |Imk_(z)L| versus the normalized frequency ωL/c for the operating mode in the 200 W, 1950 MHz, 3% bandwidth PCRB TWA, respectively. In terms of the normalized imaginary wavenumber |Imk_(z)L|, the small-signal intensity gain per axial period is expressed as

$\begin{matrix} {G = {{10\mspace{11mu}\log_{10}^{{\mathbb{e}}^{2{{{Im}\mspace{11mu} k_{z}L}}}}} \cong {8.69{{{Im}\mspace{11mu} k_{z}L}}\mspace{14mu}{\frac{dB}{Period}.}}}} & {{Eq}.\mspace{20mu} 130} \end{matrix}$

Shown in FIGS. 18 and 19 are plots of the real phase shift Rek_(z)L and the gain G versus the frequency f for the operating mode. In the 3% bandwidth from 1920 MHz to 1980 MHz, the gain is between 3.4 dB/cm and 6.6 dB/cm, which is adequate. It should be point out that the gain curve can be made flat by optimizing the RF circuit, so that the beam interacts it at a smaller value of the phase shift. This has been demonstrated in parametric design studies but will not be further discussed in this report.

As the ribbon electron beam interacts with the RF circuit, unwanted modes may be excited. Such unwanted modes could arise from the second or higher bands in the RF circuit. If not suppressed, they could cause the amplifier to self-oscillate. One promising technique to suppress unwanted modes is use of frequency-selective and mode-selective photonic crystals as described herein. There are various photonic crystals, ranging from one- to three-dimensional. For the purpose of illustration, two-dimensional dielectric square lattices are discussed.

FIG. 20 show the dispersion characteristics for 4 lowest-order TM modes in a photonic crystal consisting of a dielectric square lattice with rod radius r₀=0.9 cm, lattice constant a=4.9 cm, rod dielectric constantε₁=2.3, and background dielectric constantε₂=8.9. Labels Γ, X, and M on the horizontal axis follows the conversion in solid state physics. In this design, there is a narrow band gap of 100 MHz at 2 GHz, but no band gaps at lower or higher frequencies. This 100 MHz band gap will confine the operating mode and simultaneously allow all unwanted modes to transmit through, achieving the single-mode operation of the 200 W, 1950 MHz, 3% bandwidth PCRB TWA.

The photonic crystal design can still be optimized with larger values of dielectric constants and smaller lattice constants.

Because the PCRPB TWA is scalable to higher frequencies, wider bandwidth, and higher power output, the 1950 MHz PCRB TWA can be redesigned as a power amplifier for high-frequency (3-6 GHz) 3G wireless base stations as well as for future wireless base stations.

Although the present invention has been shown and described with respect to several preferred embodiments thereof, various changes, omissions and additions to the form and detail thereof, may be made therein, without departing from the spirit and scope of the invention. 

1. A RF amplifier comprising: a RF input section for receiving a RF input signal; a RF amplification section with at least one single-sided slow-wave structure having at least one photonic crystal; an electron ribbon beam that interacts with the RF input signal supported by said RF amplification section with at least one single-sided slow-wave structure having at least one photonic crystal so that the kinetic energy of said electron beam is transferred to the RF fields of said RF input signal, thus amplifying the RF input signal; and a RF output section that outputs said amplified RF input signal.
 2. The RF amplifier of claim 1, wherein said at least one single-sided slow-wave structure comprises metallic or dielectric rods, dots and plates.
 3. The RF amplifier of claim 1 further comprises wiggler magnets that focus said ribbon electron beam.
 4. The RF amplifier of claim 1, wherein said at least one single-sided slow-wave structure is associated with said RF interaction section.
 5. The RF amplifier of claim 1, wherein said ribbon electron beam comprises an aspect-ratio greater than unity.
 6. The RF amplifier of claim 1, wherein said at least one photonic crystal comprises one photonic crystal.
 7. The RF amplifier of claim 1, wherein said at least one photonic crystal comprises two photonic crystals.
 8. The RF amplifier of claim 2, wherein said dielectric rods comprise a two-dimensional and/or three-dimensional dielectric lattice.
 9. A method of forming a RF amplifier comprising: forming a RF input section for receiving a RF input signal; forming a RF amplification section with at least one single-sided slow-wave structure having at least one photonic crystal; and forming an electron ribbon beam that interacts with the RF input signal supported by said a RF amplification section with at least one single-sided slow-wave structure having at least one photonic crystal so that the kinetic energy of said electron beam is transferred to the RF fields of said RF input signal, thus amplifying the RF input signal; and forming a RF output section that outputs said amplified RF input signal.
 10. The method of claim 9, wherein said at least one single-sided slow-wave structure comprises metallic or dielectric rods, dots and plates.
 11. The method of claim 9 further comprises providing wiggler magnets that focus said ribbon electron beam.
 12. The method of claim 9, wherein said at least one single-sided slow-wave structure is associated with said RF interaction section.
 13. The method of claim 9, wherein said ribbon electron beam comprises an aspect ratio greater than unity.
 14. The method of claim 9, wherein said at least one photonic crystal comprises one photonic crystal.
 15. The method of claim 9, wherein said at least one photonic crystal comprises two photonic crystals.
 16. The method of claim 10, wherein said dielectric rods comprise a two-dimensional and/or three-dimensional dielectric lattice. 